Preprint typeset in JHEP style - HYPER VERSION 

FSU-HEP-050817, BNL-HET-05/24, UH-511-1079-05 



Exploring the BWCA (Bino-Wino Co-Annihilation) 
Scenario for Neutralino Dark Matter 



Howard Baer a , Tadas Krupovnickas 6 , Azar Mustafayev", Eun-Kyung Park", Stefano 
Profumo a c and Xerxes Tata d 

a Department of Physics, Florida State University Tallahassee, FL 32306, USA 
h High Energy Theory Group, Brookhaven National Laboratory, Upton, NY 11973, USA 
c Division of Physics, Mathemathics and Astronomy, California Institute of Technology, 
Mail Code 106-38, Pasadena, CA 91125, USA 

d Department of Physics and Astronomy, University of Hawaii, Honolulu, HI 96822, 
USA 



E-mail: baerQhep . f su . edu , tadasQquark . phy . bnl . gov 



mazarOhep . f su . edu,epark@hep . f su . edu, 
prof umoOcaltech . edtj . tataOphys . hawaii . edu 



Abstract: In supersymmetric models with non-universal gaugino masses, it is possible 
to have opposite-sign SU{2) and U{1) gaugino mass terms. In these models, the gaugino 
eigenstates experience little mixing so that the lightest SUSY particle remains either pure 
bino or pure wino. The neutralino relic density can only be brought into accord with the 
WMAP measured value when bino-wino co-annihilation (BWCA) acts to enhance the dark 
matter annihilation rate. We map out parameter space regions and mass spectra which are 
characteristic of the BWCA scenario. Direct and indirect dark matter detection rates are 
shown to be typically very low. At collider experiments, the BWCA scenario is typified 
by a small mass gap — ~ 20 — 80 GeV, so that tree level two body decays of 
Z2 are not allowed. However, in this case the second lightest neutralino has an enhanced 
loop decay branching fraction to photons. While the photonic neutralino decay signature 
looks difficult to extract at the Fermilab Tevatron, it should lead to distinctive events at 
the CERN LHC and at a linear e + e~ collider. 
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1. Introduction 

In i?-parity conserving supergravity models a stable neutralino {Z\) is the lightest su- 
persymmetric particle (LSP) over a large part of the parameter space of the model. A 
neutralino LSP is generally considered an excellent candidate to comprise the bulk of the 
cold dark matter (CDM) in the universe. The relic density of neutralinos in supersymmet- 
ric models can be calculated by solving the Boltzmann equation for the neutralino number 
density Q. The central part of the calculation is to evaluate the thermally averaged neu- 
tralino annihilation and co-annihilation cross section times velocity. The computation 
requires evaluating many thousands of Feynman diagrams. Several computer codes are 
now publicly |3| available to evaluate the neutralino relic density 

From its analysis of the anisotropics in the cosmic microwave background radiation, 
the WMAP collaboration has inferred that the CDM density of the universe is given by(4j, 

n CD Mh 2 = 0.113 ±0.009. (1.1) 

Since the dark matter could well be composed of several components, strictly speaking, the 
WMAP measurement only implies an upper limit on the density of any single dark matter 
candidate. Nevertheless, even this upper bound imposes a tight constraint on all models 
that contain such candidate particles, and in particular, on supersymmetric models with a 
conserved .R-parity Q. 

Many analyses have been recently performed in the context of the paradigm minimal 
supergravity model|| (mSUGRA), which is completely specified by the parameter set, 

m o ; r "i/2 ) A); tan/3 and sign(n). 

The mSUGRA model assumes that the minimal supersymmetric model (MSSM) is valid 
between the mass scales Q = Mgut and Q = M wea ^. A common value mo (m^) ((A))) 
is assumed for all scalar mass (gaugino mass) ((trilinear soft SUSY breaking)) parameters 
at Q = Mgut, and tan/3 is the ratio of vacuum expectation values of the two Higgs fields 
that give masses to the up and down type fermions. The magnitude of the superpotential 
Higgs mass term fi, but not its sign, is fixed so as to reproduce the observed Z boson mass. 
The values of couplings and other model parameters renormalized at the weak scale can be 
computed via renormalization group (RG) evolution from Q = Mgut to Q = M wea k . Once 
these weak scale parameters that are relevant to phenomenology are obtained, sparticle 
masses and mixings may be computed, and the associated relic density of neutralinos can 
be determined. 

In most of the allowed mSUGRA parameter space, the relic density Q^h 2 turns out 
to be considerably larger than the WMAP value. Consistency with WMAP thus implies 
that neutralinos should be able to annihilate very efficiently. In the mSUGRA model, the 
annihilation rate is enhanced in just the following regions of parameter space, where the 
sparticle masses and/or the neutralino composition assume special forms. 

• The bulk region occurs at low values of mo and m 1( / 2 @, §]■ In this region, neutralino 
annihilation is enhanced by ^-channel exchange of relatively light sleptons. The bulk 
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region, featured prominently in many early analyses of the relic density, has been 
squeezed from below by the LEP2 bound on the chargino mass m^ > 103.5 GeV 
and the measured value of the branching fraction B(b — > 57), and from above by the 
tight bound from WMAP. 

The stau co-annihilation region occurs at low mo for almost any m^j value where 
rrtf x ~ m z 1 - ^he s t aus > being charged, can annihilate rapidly so that f±Zi co- 
annihilation processes that maintain Z\ in thermal equilibrium with fx, serve to 
reduce the relic density of neutralinos . 

The hyperbolic branch/focus point (HB/FP) region at large mo ~ several TeV, 
where \fi\ becomes small, and neutralinos efficiently annihilate via their higgsino 



components [10]. This is the case of mixed higgsino dark matter (MHDM). 



The A-annihilation funnel occurs at large tan j3 values when 2m 7 ~ rriA and neutrali- 



nos can efficiently annihilate through the relatively broad A and H Higgs resonances [11] 
In addition, a less prominent light Higgs h annihilation corridor occurs at low m 1 / 2 |12] and 



a top squark co-annihilation region occurs at particular Aq values when m^ ~ m^; |0 

Many analyses have also been performed for gravity-mediated SUSY breaking models 
with non-universal soft terms. Non-universality of soft SUSY breaking (SSB) scalar masses 
can, 1. pull one or more scalar masses to low values so that "bulk" annihilation via i-channel 



exchange of light scalars can occur] 14, 15], 2. they can bring in new near degeneracies of 



various sparticles with the Z\ so that new co-annihilation regions open up|l6|, 15, [l7]], 3. 
bring the value of tua into accord with 2m^ so that funnel annihilation can occur [18, [0|, or 
4. they can pull the value of \x down so that higgsino annihilation can occur [18, IS, [l5|]. It is 
worth noting that these general mechanisms for increasing the neutralino annihilation rate 
can all occur in the mSUGRA model. Moreover, in all these cases the lightest neutralino 
is either bino-like, or a bino-higgsino mixture. 

If non-universal gaugino masses are allowed, then qualitatively new possibilities arise 



that are not realized in the mSUGRA model[2C, [21], 22, 23]. One case, that of mixed 
wino dark matter (MWDM), has been addressed in a previous paper [24]. In this case, 
as the weak scale value of SU{2) gaugino mass M2(weak) is lowered from its mSUGRA 
value, keeping the hypercharge gaugino mass Mi (weak) fixed, the wino component of Z\ 
continuously increases until it becomes dominant when M2(weak) < Mi (weak) (assuming 
\n\ is large). The Z1W12W coupling becomes large when Z\ becomes wino-like, resulting 
in enhanced Z\Z\ —* W + W~ annihilations. Moreover, co-annihilations with the lightest 
chargino and with the next-to-lightest neutralino help to further suppress the LSP ther- 
mal relic abundance. Indeed, if the wino component of the neutralino is too large, this 
annihilation rate is very big and the neutralino relic density falls well below the WMAP 
value. 

A qualitatively different case arises in supersymmetric models if the SSB gaugino 
masses Mi and M2 are of opposite sign. 1 As we will see below, the transition from a bino- 
like Z\ to a wino-like Z\ is much more abrupt as — M 2 (weak) passes through Mi (weak). 

The sign of the gaugino mass under RG evolution is preserved at the one loop level. 
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Table 1: Relative gaugino mass parameters at Q — Mqut and their relative values evolved to 
Q = Mz in various possible Fh irreducible representations in SU(5) and 5O(10) SUSY GUTs 
with a non-minimal GKF. Here, G442 = SU(4) x SU(2) x SU(2), and H 5 i denotes the flipped 
SU(5) x U(l) symmetry group. 



Opposite sign masses for SU(3) relative to SU(2) and U(l) gaugino mass parameters are 
well known to arise in the anomaly-mediated SUSY breaking (AMSB) model p5| . An 
opposite sign between bino and wino masses, which is of interest to us here, can arise in 
supersymmetric models with a non-minimal gauge kinetic function (GKF). In supergravity 
Grand Unified Theories (GUT), the GKF Jab must transform as the symmetric product 
of two adjoints of the GUT group. In minimal supergravity, the GKF transforms as a sin- 
glet. In SU (5) SUGRA-GUT models, it can also transform as a 24, 75 or 200 dimensional 
representation [p6| , while in 5O(10) models it can transform as 1, 54, 210 and 770 dimen- 
sional representations 1 27, S3]. Each of these non-singlet cases leads to unique predictions 
for the ratios of GUT scale gaugino masses, though of course (less predictive) combina- 
tions are also possible. The GUT scale and weak scale ratios of gaugino masses are listed 
in Table || for these non-singlet representations of the GKF. If the GKF transforms as a 
linear combination of these higher dimensional representations, then essentially arbitrary 
gaugino masses are allowed. In this report, we will adopt a phenomenological approach 
as in Ref. p4|] , and regard the three MSSM gaugino masses as independent parameters, 
with the constraint that the neutralino relic density should match the WMAP measured 
value. However, in this paper, we will mainly address the special features that arise when 
the SU{2) and U{\) gaugino masses have opposite sign. 

Much work has already been done on evaluating the relic density in models with 
gaugino mass non-universality. Prospects for direct and indirect detection of DM have also 
been studied. Griest and Roszkowski first pointed out that a wide range of relic density 
values could be obtained by abandoning gaugino mass universality p9[| . A specific form of 
gaugino mass non-universality occurs in AMSB models mentioned above, where the gaugino 
masses are proportional to the /3-functions of the corresponding low energy gauge groups: 
Mi : M2 : M3 ~ 3 : 1 : — 10. In this case the Z\ is almost a pure wino and so can annihilate 
very efficiently, resulting in a very low thermal relic density of neutralinos. This led Moroi 
and Randall |?0| to suggest that the decay of heavy moduli to wino-like neutralinos in the 
early universe could account for the observed dark matter density. Corsetti and Nath 
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investigated dark matter relic density and detection rates in models with non-minimal 
SU(5) GKF and also in O-II string models[31]. Birkedal-Hanson and Nelson showed that 
a GUT scale ratio M1/M2 ~ 1.5 would bring the relic density into accord with the measured 
CDM density via MWDM, and also presented direct detection rates [|3^]. Bertin, Nezri and 
Orloff studied the variation of relic density and the enhancements in direct and indirect DM 
detection rates as non- universal gaugino masses are varied 1 33]. Bottino et al. performed 
scans over independent weak scale parameters to show variation in indirect DM detection 
rates, and noted that neutralinos as low as 6 GeV are allowed [34]. Belanger et al. have 
recently presented relic density plots in the too vs. toi/ 2 plane for a variety of universal and 
non-universal gaugino mass scenarios, and showed that large swaths of parameter space 
open up when the SU(3) gaugino mass M3 becomes small |j3"5| : this is primarily because 
the value of |/i| reduces with the corresponding value of M3, resulting in an increased 
higgsino content of the neutralino. Mambrini and Muhoz, and also Cerdeno and Muhoz, 
examined direct and indirect detection rates for models with scalar and gaugino mass non- 
universality |3^]. Auto et a/.|i~6[| proposed non-universal gaugino masses to reconcile the 
predicted relic density in models with Yukawa coupling unification with the WMAP result. 
Masiero, Profumo and Ullio exhibit the relic density and direct and indirect detection 
rates in split supersymmetry where Mi, M2 and fi are taken as independent weak scale 
parameters with ultra-heavy squarks and sleptons ]j37|l. 

The main purpose of this paper is to examine the phenomenology of SUSY models 
with non-universal gaugino masses and examine their impact upon the cosmological relic 
density of DM and its prospects for detection in direct and indirect detection experiments, 
and finally for direct detection of sparticles at the Fermilab Tevatron, the CERN LHC and 
at the future international e + e~ linear collider (ILC). Towards this end, we will adopt a 
model with GUT scale parameters including universal scalar masses, but with independent 
SU(2) and U(l) gaugino masses which can be of opposite-sign. Indeed while some of the 
earlier studies with non-universal gaugino masses mentioned in the previous paragraph do 
allow for negative values of Mi /Ma, we are not aware of a systematic exploration of this 
part of parameter space. For the most part, we adjust the gaugino masses until the relic 
density matches the central value determined by WMAP. Whereas the case of same-sign 
gaugino masses allows consistency with WMAP via both bino-wino mixing and bino-wino 
co-annihilations (the MWDM scenario) , the opposite sign case admits essentially no mixing 
between the bino and wino gaugino components. Agreement with the WMAP value can 
be attained if the LSP is bino-like, and the wino mass M2 — —Mi at the weak scale, so 
that bino-wino co-annihilation (BWCA) processes come into play and act to reduce the 
bino relic density to acceptable values. 

The BWCA scenario leads to a number of distinct phenomenological consequences. 
For direct and indirect DM search experiments (except when sfermions are also very light), 
very low detection rates are expected in the BWCA scenario because SU {2)l x U(1)y gauge 
invariance precludes couplings of the bino to gauge bosons. Regarding collider searches, 
the BWCA scenario yields relatively low Z2 — Z\ mass gaps, so that two-body tree level 
neutralino decays are not kinematically allowed. However, when too ~ .5 — 1 TeV, the 
loop-induced radiative decay Z<i — > Zij is enhanced, and can even be the dominant Z2 
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decay mode. This may give rise to unique signatures involving isolated photon plus jet(s) 
plus lepton(s) plus E™ ms events at the Fermilab Tevatron and the CERN LHC hadron 
colliders. At the ILC, Z1Z2 production can lead to "photon plus nothing" events, while 
Z2Z2 production can lead to diphoton plus missing energy events at large rates. It would 
be interesting to examine if these signals can be separated from SM backgrounds involving 
neutrinos, e + e~7 and multiple gamma processes where the electron and positron, or some of 
the photons, are lost down the beam pipe. Potentially, the energy spectrum of the signal 
7 may allow for the extraction of and via photon energy spectrum endpoint 

measurements. 

The remainder of this paper is organized as follows. In Sec. ^, we present the param- 
eter space for the BWCA scenario, and show the spectrum of sparticle masses which are 
expected to occur. We also discuss some fine points of the BWCA relic density analysis. 
In Sec. |3], we show rates for direct and indirect detection of DM in the BWCA scenario. 
These rates are expected to be below detectable levels unless some sfermions are very light 
or additional annihilation mechanisms can be active. In Sec. |3], we present expectations 
for the radiative neutralino decay Z<i — ► Z\~y in BWCA parameter space. In Sec. [5], we 
examine the implications of the BWCA scenario for the Fermilab Tevatron, CERN LHC 
and the ILC. In Sec. ^, we present our conclusions. In the Appendix we adapt the idea 
of integrating out heavy degrees of freedom, familiar in quantum field theory, to quantum 
mechanics, and use the results to obtain simplified expressions for neutralino masses in the 
large limit where the higgsinos can be "integrated out". 

2. Sparticle mass spectrum in the BWCA scenario 

It is well known that a pure bino LSP can annihilate rapidly enough to give the observed 
relic density only if scalars are sufficiently light, as for instance in the so-called bulk region 
of the mSUGRA model. Our goal here is to explore SUGRA models with universal values of 
high scale SSB scalar masses and ^4-parameters, but without the assumption of universality 
on gaugino masses. In models without gaugino mass universality, the annihilation rate of 
a bino LSP may be increased in several ways, including 

1. by increasing the higgsino content of the LSP, which may be achieved by decreasing 
the gluino mass relative to the electroweak gaugino masses J 



2. by increasing the wino content of the LSP, by reducing the ratio M1/M2 as in the 



MWDM scenario which has been the subject of many investigations |32|, 24] 



3. by allowing co-annihilations between highly pure bino-like and wino-like states with 
comparable physical masses. We dub this the bino-wino coannihilation (BWCA) 
scenario. 

The MWDM scenario is realized when the SU(2) gaugino mass parameter M^weak) = 
M2(M wea k) approaches the U(l) gaugino mass parameter Mi (weak). Since the Z\ in 
MWDM can annihilate very efficiently into W + W~ pairs via its wino component, care 
must be taken to ensure that this wino component is not so large that the relic density 
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falls below its WMAP value. 2 The reader can easily check that if the electroweak gaugino 
masses are equal at the weak scale, the photino state 



7 = cos 9wB + sin OyyW, 



is an exact mass eigenstate of the tree level neutralino mass matrix, and has a mass equal 



to the common weak scale gaugino mass, M [38]. Since the photino can readily annihilate 
to W + W~ pairs via chargino exchange in the t channel, we would expect that it gives a 
relic density that is smaller than the WMAP value: in other words, the wino content of the 
Z\ must be smaller than sin Oy/ in order to obtain the WMAP value of the relic density. 
We will see shortly that this is indeed the case. 

In the MWDM scenario, treating the gaugino-higgsino mixing entries in the neutralino 
mass matrix, whose scale is set by Myy, as a perturbation, the signed tree level gaugino 
masses for the case Mi (weak) = M2(weak) = M are given (to second order) by, 



M + 



M (photino) 
(sin/3 — cos/3) 2 (sin/3 + cos/3) 2 



M 



+ 



M + fj, 



(zino) 



(2.1) 



If instead of exact degeneracy between Mi (weak) and M2(weak), we have Mi (weak) = 
M — 5 and M2(weak) = M + 5, where 5 is as small or comparable to the gaugino-higgsino 
mixing entries, these eigenvalues change to, 



where 



and 



1 



Jpcrt 



M 2 



2 g 2 + g 



12 



(sin /3 — cos /3) 2 (sin (5 + cos /3) 



M 



/I 



+ 



M + (JL 



2 - {g 2 + g' 2 ) 2 + % 2 - ^ 



+ 4- 



5 2 



(2.2) 



(2.3) 



(2.4) 



■"pert v peit 

Clearly, these eigenvalues reduce to the masses (^J) of the photino and zino states when 
5^0. 

The generic case where the differences between Mi (weak), M2(weak) and ±// are all 
much larger than the gaugino-higgsino mixing entries is much simpler to treat since we do 
not have to worry about degeneracies as in the Mi (weak) ~ M2(weak) case that we have 
just discussed. In the absence of gaugino-higgsino mixing, the bino and neutral wino are the 
gaugino mass eigenstates. Then, again treating the gaugino-higgsino mixing entries in the 
neutralino mass matrix as a perturbation, we see that mixing between wino and bino states 
occurs only at second order in M\y/M, where M denotes a generic mass difference between 
the "unperturbed" higgsino and gaugino mass eigenvalues. This is in sharp contrast to the 
MWDM scenario where even the "unperturbed" gaugino states are strongly mixed because 



2 We emphasize that, while a neutralino relic density smaller than the WMAP value is not excluded, in 
this paper we confine ourselves to scenarios that accommodate this value. 
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of the degeneracy of the eigenvalues when Mi (weak) = M2(weak). It is simple to show 
that the signed masses of the bino-like and wino-like neutralino eigenstates take the form, 



Mi + -M| sin 2 6 W 
M 2 + -M| cos 2 9 W 



(sin — cos (3) 2 ^ (sin f3 + cos /3) 5 



Mi - n Mi+ n 

(sin/3 — cos/3) 2 (sin/3 + cos/3) 21 
M 2 -n ' M 2 + m 



(bino), 

(wino), (2.5) 



where Mi and M 2 in ( |2.5| ) denote their values at the weak scale. Note that for the special 
case Mi (weak) = — M2(weak), the formulae Q2.5| ) apply even though the physical masses 
of the two lighter neutralinos may be very close to one another. The main purpose of 
our discussion of the approximate masses and mixing patterns is that these will enable 
us to better understand the differences in phenomenology in the MWDM scenario, where 
Mi (weak) is slightly smaller than M2(weak), and the BWCA scenario where Mi (weak) ~ 
-M 2 (weak). 

Before proceeding to explore other sparticle masses in these scenarios, we first check 
that the qualitative expectations for neutralino mixing patterns and relic density discussed 
above are indeed realized by explicit calculation. Toward this end, we adopt the subpro- 
gram Isasugra, which is a part of the Isajet v 7.72' event generator program|3^]. 3 Isasugra 
allows the user to obtain sparticle masses for a wide variety of GUT scale non-universal 
soft SUSY breaking terms. The sparticle mass spectrum is generated using 2-loop MSSM 
renormalization group equations (RGE) for the evolution of all couplings and SSB parame- 
ters. An iterative approach is used to evaluate the supersymmetric spectrum. Electroweak 
symmetry is broken radiatively, so that the magnitude, but not the sign, of the superpoten- 
tial [i parameter is determined. The RG-improved 1-loop effective potential is minimized 
at an optimized scale to account for the most important 2-loop effects. Full 1-loop radia- 
tive corrections are incorporated for all sparticle masses. To evaluate the neutralino relic 



density, we adopt the IsaRED program|3[, which is based on CompHEP[40| to compute 
the several thousands of neutralino annihilation and co-annihilation Feynman diagrams. 
Relativistic thermal averaging of the cross section times velocity is performed The 
parameter space we consider is given by 

"^0, ?7ii/2) Ab tan/3, sign(fi), M\ or M2, (2-6) 

where we take either Mi or M2 to be free parameters (renormalized at Q = Mgut), and 
in general not equal to ^1/2 • When Mi is free, we maintain M2 = Tnx/ 2 , while when M2 is 
free, we maintain Mi = m-i/2- 

We show our first results in Fig. [l], where we take tuq = = 300 GeV, with Aq = 0, 
tan/3 = 10, /J, > with rat = 178 GeV. We plot the neutralino relic density ^j^ 2 in frame 
a) versus variation in the U(l) gaugino mass parameter Mi. For Mi = 300 GeV, we are 
in the mSUGRA case, and £1% h 2 = 1-3, so that this model would be strongly excluded 
by the WMAP measurement. For smaller values of Mi, the bino-like neutralino becomes 



3 Isajet 7.72' is Isajet v7.72 with the subroutine SSM1LP modified to yield chargino and neutralino mass 
diagonalization all at a common scale Q = v /mg 1 • 
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lighter and two dips occur in the neutralino relic density. These correspond to the cases 
where 2m^ — m-h and Mz as one moves towards decreasing Mi, i.e. one has either light 
Higgs h or Z resonance annihilation 4 . Instead, as M\ increases past its mSUGRA value, 
the Z\ becomes increasing wino-like and the relic density agrees with the WMAP value at 
Mi ~ 490 GeV, where the MWDM scenario is realized [24]. For yet larger values of Mi, the 
wino content of Z\ is so large that the relic density falls below the WMAP value. Turning 
to negative values of Mi, we see that as Mi starts from zero and becomes increasingly 
negative, the Z and h poles are again encountered. The relic density is again much larger 
than the WMAP bound until it begins decreasing for Ml < -400 GeV. At Mi ~ -480 
GeV, the relic density is again in accord with the WMAP value, while for more negative 
values of Mi, the relic density is too low, so that some other form of CDM or non-thermal 
production of neutralinos would be needed to account for the WMAP measurement. 

In frame b), we show the amplitude, i2g^, for the bino/wino content of the lightest 

neutralino Z\. Here, we adopt the notation of Ref. wherein the lightest neutralino 

is written in terms of its (four component Majorana) Higgsino and gaugino components as 

Zi = v ( i ] i> h ° + + 4 1} A 3 + vjpXo, (2.7) 

where R^ = \v^\ and = |. A striking difference between the positive and negative 
Mi portions of frame b) is the shape of the level crossings at R^ = R^ = 1/ y/2: while the 
transition from a bino-like to a wino-like LSP is gradual when Ml > 0, it is much more 
abrupt for negative values of Mi. We had already anticipated this when we noted that 



4 For values of \M X \ ~ 100 GeV, we have checked that the contribution of F(Z — > ZiZ^) is always below 
limits from LEP on the invisible width of the Z boson. 




Figure 1: A plot of a) relic density flcDMh 2 and 6) bino/wino component of the lightest neutralino 
as a function of M\ for mo = 300 GeV, mi/2 = 300 GeV, Aq = 0, tan f3 = 10, fi > and mt = 178 
GeV. 
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for Mi (weak) = M2(weak) the mass eigenstates are the photino, and aside from a small 
admixture of higgsinos, the zino, while the corresponding eigenstates were bino-like and 
wino-like as long as | M\ (weak) — M2(weak)| was larger than the gaugino-higgsino mixing 
entries in the neutralino mass matrix. We also see that, for M\ > 0, the MWDM scenario 
is realized for the value M\ where < sm8\y — 0.48 as we also expected. Finally, again 
as we anticipated, for Mi < 0, the BWCA scenario is obtained for |Mi| just above the level 
crossing, when the LSP is mainly bino-like and close in mass to the wino-like Z\. For a 
bino-like LSP, annihilation to vector boson pairs is suppressed, and the only way to reduce 
the relic density in the case of large negative Mi is by having rather small Z 2 — Z\ and 
W\ — Z\ mass gaps, so that chargino and neutralino co-annihilation effects are large |p3| , 
and act to decrease the relic density. Of course, when —Mi > M2, then the Z\ suddenly 
becomes pure wino-like, and very low relic density is obtained, as in the case of AMSB 
models. 

Similar results are obtained by keeping Mi = m 1( / 2 , and varying M2, as shown in Fig. 
|2[ The MWDM scenario is reached just to the right of the level crossing where M2(weak) is 
slightly larger than Mi (weak), while the BWCA scenario is reached for M2 = —156 GeV, 
where |M2(weak)| is again just above | Mi (weak) | and the Z\ is essentially a bino. Again, 
the Z 2 — Z\ and, as discussed just below, also the W\ — Z\, mass gaps are very small, so 
that co-annihilation plays an essential role in reducing Qg^h 2 to an acceptable value. We 
see also that there is a region of mainly small negative M2 values where < , where 
a charged LSP is obtained. This is, of course, ruled out by the negative results for searches 
for charged stable relics from the Big Bang. 

In both the MWDM and the BWCA scenarios, Z\, Z2 and W\ are the lightest spar- 
ticles. They will likely have a large impact upon the phenomenology of these models. Of 
particular interest is the ordering of the mass spectrum of these particles. We work this out 
in the so-called large |/j| approximation, where |//| 3> Mi,M2 applicable in many models. 
The chargino sector, since it consists of just two states, is simple and is given by a 
relatively simple and well known expression that we do not reproduce here (see e.g. Ref. 
]3q1 ). The neutralino sector is much more complicated but, in the large approxima- 
tion, it is possible to "integrate out" the higgsinos, and work with an effective low scale 
theory that only includes the neutral bino and wino as discussed in the Appendix. The 
couplings of winos and binos to higgsinos in the original theory manifests itself as a mix- 
ing between winos and binos in the effective theory, where this mixing is suppressed by 
1/|/j|. To 0(l//i 2 ), the signed tree level neutralino masses are given (in terms of weak scale 
parameters) by, 5 

. M^ M 1 + M 2 + asec 2 9 w + M 2 r ] ±K 1 

V /i 2 7 Z L2 1(1+??) 

where 

Mh sin 2(3 . „ . 
a = — ^ < Mi , M 2 

A* 

We are abusing notation here in that we are using and m^ 2 to denote the signed neutralino mass, 
whereas everywhere else we use the same symbols to denote the physical (positive) neutralino masses. 
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V= -l + tan 2 % < 1, 

and 

K 2 = (Mi - M 2 ) 2 + 2a (Mi + M 2 ) sec 2 (% + a 2 sec 4 H - 
+ 2 (Mi + M 2 )M 2 r] - 4a(Mi + M 2 tan 2 W ) - 4MiM 2 77- 



We emphasize that ( |2.8[ ) is valid even when the gaugino masses are comparable to the 
gaugino-higgsino mixing terms in the neutralino mass matrix, and is a useful approximation 
as long as \fi\ is large. If Mi is not very close to M 2 , we can simplify the expressions for 
m Zi 2 ^ expanding k ~ \M\ — M 2 | plus terms that are suppressed by powers of fi. Such 
an expansion then yields, 

M w 2n . nn , r M 2 ^ sin 2 2/?tan 2 6liy . o N 

m~ = Mi + ^tan 2 ^sm2/3-Mi^ / *, + ° (Vfl ). 

^ Mi — M 2 

™z 2 = M 2 + ^ sin 2/3 - M 2 ^f _ sin 2 2/3 + 0(1/^), (2.9) 

2 fi fi z fi z Mi — M 2 

M 2 M 2 
= M 2 + — ^ sin 2/3 - -^M 2 + 0(1/ fi 3 ) , 

where we have assumed \M%\ < |M 2 | when we associated the bino-like state with Z\. 
Except for this, (|2.9|) are valid for all magnitudes and signs of Mi and M 2 , as long as 



|Mi — M 2 | ^> Myy/\fi\- In particular, we can use these expressions for the BWCA scenario, 
but for the MWDM case, we would have to use fl2.g| ) to get the neutralino masses. 

The ordering of the mass spectrum of the lightest SUSY particles in the same sign 
MWDM and opposite sign BWCA scenarios is different, as can be seen from Fig. 0, where 



m =300 GeV, m 1/2 =300 GeV, tanp=10, A =0, u>0, m=178 GeV 
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Figure 2: A plot of a) relic density flcDMh 2 and 6) bino/wino component of the lightest neutralino 
as a function of Mi for mo = 300 GeV, m\/2 — 300 GeV, Aq = 0, tan j3 = 10, fi > and mt = 178 
GeV. 
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Figure 3: The mass splitting between the lightest chargino and the next-to-lightest neutralino, 
m w 1 ~ m z 2 ^ n BWCA scenario (dot-dashed blue line, with opposite M\ and signs) and in the 
MWDM scenario (dashed red line, with same M\ and M2 signs), as a function of the absolute GUT 
scale value of |Mi|/mi/2- The arrows indicate the value of \M\\/m\/2 where the neutralino thermal 
relic abundance equals the central WMAP CDM abundance value. The other SUSY parameters 
are as in Fig. 1. 

we plot the W1-Z2 mass splitting as a function of the absolute GUT-scale value \Mx\/m\/2- 
For the MWDM case, the lightest chargino is lighter than the next-to-lightest neutralino, 
while the opposite holds true in the BWCA case. As long as | M\ (weak) | < |M2(weak)|, 
both W\ and Zi are dominantly wino-like, and the tree-level mass splitting between them 
can be read off from (|2.9| ), as long as the weak scale values of M\ and M2 are not too close. 
We then find that at tree- level 44], 



ml^, tan 9\y , 

Am tree = m^-. — rriy & rr tit-, \ s ox sin 2p. (2.10) 

w ^ z * (Mi (weak) - M 2 (weak) n 2 ) 



for all combinations of signs of gaugino masses. For the BWCA case, the denominator 
I Mi (weak) — M2(weak)| is very large, so that the tree level splitting is negligible compared 
to the one-loop splitting; the latter is always positive, hence in the opposite sign BWCA 
case, when the LSP is bino-like, the lightest chargino is always heavier than the next-to- 
lightest neutralino. In contrast, for the same-sign MWDM case, the tree level splitting 



(over the range of Mi values where (2.1C) is valid) is negative, and comparable to or larger 
than the one-loop splitting, so that the chargino is now usually heavier than Z2, as may 
be seen by the dashed line in Fig. ||. 

Various other sparticle masses are also affected by varying the gaugino masses, since 
these feed into the soft term evolution via the RGEs. In Fig. El we show the variation of 
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the sparticle mass spectrum with respect to the GUT scale ratio M\jmxj 2 for the same 
parameters as in Fig. |]. In the mSUGRA case where M\/m,\/2 = 1, there is a relatively 
large mass gap between Z 2 and Z\\ m z 2 ~ m z 1 = 106.7 GeV. As M\ varies to large 
positive values (the MWDM case), the mass gap shrinks to ra^ - = 31.9 GeV. As 
Mi varies to large negative values, the mass gap also decreases, this time to just 22.7 GeV 
in the BWCA scenario. We also note that as \M\\ increases, the Br, /zr and f\ masses 
also increase, since M\ feeds into their mass evolution via RGEs. This also gives rise to 
the nearly symmetric behavior versus the sign of M\ for the mass spectrum of first and 
second generation sfermions. As the coefficient appearing in front of M\ in the RGEs is 
larger (and with the same sign) for the right handed sfermions than for the left handed 
ones, one expects, in general, a departure from the usual mSUGRA situation where the 
lightest sleptons are right-handed. As a matter of fact, whereas in mSUGRA mg L >> mi R 
for rriQ ~ mi/2) i n the case of BWCA or MWDM, we find that mg L ~ m^ R . As shown 
in the figure, the right-handed squark masses also increase with increasing \M±\, although 
the relative effect is less dramatic than the case involving sleptons: the dominant driving 
term in the RGEs is, in this case, given by M3 (absent in the case of sleptons), hence 
variations in the GUT value of M\ produce milder effects. The trilinear SSB At, Ab and 
A T parameters have a linear dependence on gaugino mass in their RGEs, which means the 
weak scale ^4-parameters will be asymmetric versus the sign of M\. The \i parameter is 
also slightly asymmetric. This gives rise to the asymmetric behavior of the third generaton 
sfermion masses with respect to the sign of the gaugino mass. 

In Fig. H we show a plot of sparticle masses for the same parameters as in Fig. |], but 
versus M^/min^ In this case, as \M 2 \ is decreased from its mSUGRA value of 300 GeV, 

the W\ and Z 2 masses decrease until h 2 reaches 0.11 in both the BWCA and MWDM 

^1 

scenarios. In this case, with decreasing I-M2I, the left- slepton and sneutrino masses also 
decrease, again leading to mg L ~ mg R . The left-handed squark masses similarly decrease. 
This increase is more pronounced that in Fig. || because the SU{2) gauge coupling is larger 
than the U(l)y gauge coupling. The SU(2) singlet right-handed sfermion masses are not 
affected, with the net result that the mSUGRA mg L >> rrie R hierarchy is again altered. 

It should be apparent now that most points in the niQ vs. mj/j plane can become 
WMAP allowed by adopting an appropriate negative value of either Mi or M2 such that 
one enters into the BWCA scenario. The exception occurs if the WMAP-allowed point 
is obtained because the ^4-funnel or stau co-annihilation region is reached instead. To 
illustrate this, we plot in Fig. ^ the ratio r\ = M\/mi/ 2 in frame a) or r 2 = M 2 /mi/ 2 in 
frame b) needed to achieve a relic density in accord with the WMAP central value. We 
see in frame a) that r\ generally increases as one moves from lower-left to upper-right. 
The structure in the upper-left of the plot occurs when —Mi is dialed to such a value that 
2m^ ~ niA, i-e. one is entering the A-funnel (even though tan/3 is relatively low) instead 
of the BWCA scenario. These regions will of course have a much larger Z 2 — Z\ mass 
gap than points in the BWCA scenario. Like the non-universal mass scenario [15|, the 
BWCA scenario allows the A funnel to be reached for any value of tan f3, but should be 
distinguishable from this because the Z\ — Z 2 mass gap, for instance, will be quite different 



- 12 - 



m g =300GeV, m 1/2 =300GeV, tan p =10, A g =0, [i >0, m, =178GeV 
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Figure 4: A plot of various sparticle masses vs. Mi/m 1 / 2 for mo = 300 GeV, m 1 / 2 = 300 GeV, 
A = 0, tan/3 = 10 and /j > 0. 



in the two scenarios. In frame b), the ratio r% that gives rise to a WMAP-allowed point 
is shown to increase as one travels from lower to higher values of rrii/2- In this case, since 
the value of M<i hardly changes the value of ni^^, the A- funnel is never reached, and the 
BWCA region can be accessed over most of parameter space, save near the left-hand edge 
in the stau co-annihilation region. 

3. Dark matter in the BWCA scenario 

3.1 Neutralino relic density in the BWCA scenario: a closer look 

In order to better understand the coannihilation mechanisms which drive the neutralino 
relic abundance within the WMAP preferred range in the BWCA vs. the MWDM scenario 
(with opposite and same signs for Mi and M2), we adopt the sample point defined by the 
mSUGRA input parameters (mo=300 GeV, m 1 / 2 =300 GeV, tan/3=10, sgn(/i) > 0, A)=0, 
"i-top=178 GeV), and pick the Mi/m 1 / 2 values which give = 0.11, i.e., respectively, 

Mi/mi/2 = —1.619 and Mi/m 1/ / 2 = 1.630. We plot in Fig. the thermally averaged cross 
section including coannihilations (solid) and without coannihilations (dashed lines) times 
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m g =300GeV, m 1/2 =300GeV, tan$=10, A g =0, n >0, m t =178GeV 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

M 2 lm v2 

Figure 5: A plot of various sparticle masses vs. M2/m 1 / 2 for mo = 300 GeV, m 1 / 2 = 300 GeV, 
A = 0, tan/3 = 10 and \l > 0. 

the relative velocity as a function of the temperature T. In the same sign MWDM case, 
one notices that the very significant wino-component in the photino-like lightest neutralino 
gives a non- negligible s-wave contribution via annihilation to W pairs, (i.e., a contribution 
which goes as (crv)(T) ~ a), while in the opposite sign BWCA case the lightest neutralino 
is a pure bino, and the squark mediated i-channel dominated pair annihilation cross section 
is strongly s-wave suppressed ({crv)(T) ~ a' + b^T/m^sp), a' <C b'). As a result, the pair 
annihilation cross section at T = (relevant for indirect DM detection) is more than one 
order of magnitude suppressed in the BWCA case. 

A second effect, also traced back to the lack of a significant wino component in the 
opposite sign BWCA case, is the role and onset of coannihilations. To isolate the role of 
coannihilations, we show by dashed lines (av)(T) computed without any coannihilation 
contribution. In the BWCA case, coannihilations play a much more important role, as can 
be understood by looking at the relative size of the cross sections with and without coanni- 
hilations around the lightest neutralino freeze-out temperature Xf. ., indicated by an arrow 
in the figure. Furthermore, in the BWCA case, the onset of the coannihilation regime takes 
place at lower temperatures, a reflection of the reduced mass splitting between the coan- 
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nihilation partners and the LSP. At the freeze-out temperature, {av){T) is dominated by 
coannihilations in the opposite sign BWCA case, while the relative coannihilation contribu- 
tion for the MWDM case (for which the Z\s can annihilate to WW) is significantly smaller. 
Observe also the indicated relic abundance of the LSP computed without coannihilations, 
respectively 0.759 in the BWCA case and 0.138 in the same MWDM. 

Coannihilations of particles A and B effectively enter the neutralino pair annihilation 
cross section when the center of mass momentum p c .m. of the neutralino-neutralino system 
satisfies the relation 



s = Ap 2 cra + 4?™^ > (m A + m B ) 2 - 



(3.1) 



A suitable quantity to illustrate the onset of coannihilations is given by an effective anni- 
hilation rate (av) e fi, defined as in Ref. [45| 



(<™) e ff 



W cS (p c . m .) 



+ m ~ 



(3.2) 



(we refer the reader to Ref. [p| for the definition of the effective annihilation rate W e s, 
which is essentially the sum over all (co-) annihilation channels, properly weighted, of the 
various annihilation rates per unit volume and unit time) and such that 



lim (av) eS = (av)(T = 0). 

Pc.m.— »0 



(3.3) 



BWCA: tan$=10, A =0, n >0, m, =178 GeV 



BWCA: M 2 * m 1/2 , tanf>=10, A =0, |i >0, m, =178 GeV 




Figure 6: Contours of a) r\ and b) V2 in the mo vs. mi/ 2 plane for tan/3 = 10, Aq = 0, /i > 0. 
Each point has h 2 — 0.11. The origin of the strips in the LEP excluded region of the right hand 
frame is discussed in Sec. pi 
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Figure 7: The thermally averaged effective pair annihilation cross section times velocity, (av) (T) 
as a function of the temperature, for the BWCA (black lines) and MWDM (red lines) scenarios. 
The dotted lines correspond to (av){T) without the inclusion of coannihilation effects. The arrow 
indicates the approximate temperature where freeze-out occurs. 



The temperature dependence of (av)(T) is factored out in the weight function K(p c , m ,,T) 



1 45], so that 



M(T) = fdfa. ^N «( P , ni .,T). (3.4) 

In Fig. [8|, we plot the two effective annihilation rates (crv) c f[ for the two BWCA and 
MWDM sample cases at the mSUGRA point given above, with the correct WMAP relic 
abundance. The coannihilation thresholds are shifted to larger p c . m . values in the same sign 
case (the mass splitting of coannihilating partners is increased). The W1-Z2 mass splitting 
is too small to resolve the separate contributions, and the two bumps correspond to the 
onset of Z\-Wi, Z2 coannihilations and the onset of (co-)annihilations amongst W\ and Z2. 
When the bino-wino mixing is suppressed, the second processes contribute much more than 
the first processes to {av) c &. However, as shown by the weight function K(p c . m .,T = Tf. .), 
(av) e ff is largely sampled in a p c . m . range where the W\,Z2 (co-) annihilations are not 
kinematically accessible. This picture, however, depends on the LSP mass : had we picked 
a larger value for m^, the mass splitting between Z\ and Z2, W2 needed to obtain a 
sufficiently low relic abundance would have been smaller, so that the coannihilation bumps 
in W e ff/AE 2 would have occurred at smaller center-of-mass momenta, where the sampling 
function k has not yet fallen to very small values. In this case, the role of W\ — Z2 




Figure 8: The effective (av) e g term, as defined in Eq. (|3.2| ), as a function of the center of mass 
momentum p e g, for the BWCA (black lines) and MWDM (red lines) scenarios, with (solid lines) 
and without (dotted lines) coannihilations. The green dotted line indicates the "weight function" 
K(p c . m .,T = Tf. .). 

coannihilations would have been enhanced even in the MWDM case. 

In the limit in which sfermions are heavy, and thus the bino annihilation cross section 
is extremely suppressed, and in which \i S> Mi (weak), M2(weak) (in the remainder of 
this paragraph, we will implicitly use = Mi^weak)), i.e., in the pure gaugino limit, 
the relic abundance should only depend on i.) the LSP mass scale and on ii.) the LSP- 
wino system splitting relative to the LSP mass. Since a pure wino-like system has a relic 



abundance which goes like [46] 



/ Mo V 

ftpurc winos/l 2 ^ C • ( J , C ~ 0.024, 7 ~ 1.9, (3.5) 

the relic abundance of a coannihilating bino will be given by that of the pure wino system, 
rescaled by the exponential factor exp[— ((M2 — M\)/M\) ■ x)], with x ~ Mi/Tf . Q ., and 
rescaled by the new parasite bino degrees of freedom 6 . The wino system carries 2+4 degrees 
of freedom, while the bino 2, hence we expect an enhancement of the relic abundance for a 
coannihilating bino of a factor (4/3) 2 Ww - The relic abundance should then take the form 

~ M 2 - Mi 



"bino+wino/i ~ a ■ \ - T ^ ) • exp 



Mi 



c(4/3) 2 . (3.6) 



6 The degrees of freedom should be as well weighted according to the mass splitting, but this is a higher 
order effect. 
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Figure 9: Iso-level curves at given values of the neutralino relic abundance in the M\, (Mi 
|Mi|)/|Mi| plane. The green strip indicates the WMAP 2-a range. We assumed here (i,m^,mji 
100 • M2, trig being all sfermion masses. 



Since M2 — Mi, this allows in principle to define a strip in the Mi, M2 ^^ plane of WMAP 
preferred relic abundance. We show the result in Fig. ^, taking [A,mg,mA = 100 • M%, mg 
being all sfermion masses. Indeed, the iso-level curves for the relic abundance show the 
functional form of Eq. ( |3.6[) . 

3.2 Direct and indirect detection of neutralino CDM 

In this section, we turn to consequences of the BWCA scenario for direct and indirect 
detection of neutralino dark matter p8|. We adopt the DarkSUSY code^jj, interfaced to 
Isajet, for the computation of the various rates, and resort to the Adiabatically Contracted 
N03 Halo model[|50| for the dark matter distribution in the Milky Way 7 . We evaluate the 
following neutralino DM detection rates: 



Direct neutralino detection via underground cryogenic detectors [54]. Here, we com- 
pute the spin independent neutralino-proton scattering cross section, and compare 
it to expected sensitivities [^5| for Stage 2 detectors (CDMS2|Q, Edelweiss2||57|, 
CRESST2||, ZEPLIN2||) and for Stage 3, ton-size detectors (XENON Q, GERDA0, 
ZEPLIN4[^] and WARP (6^]). We take here as benchmark experimental reaches of 
Stage 2 and Stage 3 detectors the projected sensitivities of, respectively, CDMS2 and 
XENON 1-ton at the corresponding neutralino mass. 



7 For a comparison of the implications of different halo model choices for indirect DM detection rates, 
see e.g. Refs. ||, || [FJ|. 
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Indirect detection of neutralinos via neutralino annihilation to neutrinos in the core of 



the Sun[S4|. Here, we present rates for detection of — ► fj, conversions at Antares[55] 
or IceCube|36|]. The reference experimental sensitivity we use is that of IceCube, with 
a muon energy threshold of 25 GeV, corresponding to a flux of about 40 muons per 
km 2 per year. 

• Indirect detection of neutralinos via neutralino annihilations in the galactic center 
leading to gamma rays|37]], as searched for by EGRET f58|, and in the future by 
GLAST[69[. We evaluate the integrated continuum 7 ray flux above a = 1 GeV 



threshold, and assume a GLAST sensitivity of 1.0x10 10 cm 2 s . 

• Indirect detection of neutralinos via neutralino annihilations in the galactic halo 
leading to cosmic antiparticles, including positrons |70| (HEAT|7l"[|, Pamela ]72[ and 
AMS-02|7|), antiprotons|4| (BESS(7|, Pamela, AMS-02) and anti-deuterons (Ds) 
(BESS|76|], AMS-02, GAPS|77]]). For positrons and antiprotons we evaluate the aver- 
aged differential antiparticle flux in a projected energy bin centered at a kinetic energy 
of 20 GeV, where we expect an optimal statistics and signal-to-background ratio at 
space-borne antiparticle detectors |53|, fT8|. We take the experimental sensitivity that 
of the Pamela experiment after three years of data-taking as our benchmark. Finally, 
the average differential antideuteron flux has been computed in the 0.1 < Tq < 0.4 
GeV range, where Tq stands for the antideuteron kinetic energy per nucleon, and 



compared to the estimated GAPS sensitivity |77|| (see Ref. [79] for an updated discus- 



sion of the role of antideuteron searches in DM indirect detection) 



In Fig. 1C, we show various direct and indirect DM detection rates for mo = = 300 
GeV, with Aq = 0, tan (3 = 10 and ^ > 0, while Mi is allowed to vary. The Mi value 
corresponding to the mSUGRA model is denoted by a dot-dashed vertical line, while the 
BWCA and MWDM scenarios with Sl^h 2 = 0.11 are denoted by dash-dash-dot and dashed 
vertical lines, respectively. The dotted lines correspond to the sensitivity level of each of 
these experiments; i.e., the signal is observable only when the model prediction is higher 
than the corresponding dotted line. While the minimum sensitivity for the direct detection 
rates in frames b) - f) refers to the minimum magnitude of the signal that is detectable 
(and hence independent of the LSP mass), the smallest detectable cross section shown by 
the dotted curves in frame a) depends on the value of . 

In frame a), we plot the spin-independent neutralino-proton scattering cross section. 
We see that as Mi is decreased, and becomes increasingly negative, the neutralino-proton 
scattering cross section plummets to values in the 10 -12 pb range, far below the sensitivity 
of any planned detector. The drop-off is due to increasing negative interference amongst 
the contributing Feynman diagrams. 

In frame b) , we show the flux of muons from neutralino pair annihilations in the core of 
the Sun. The muon flux is below the reach of IceCube in the mSUGRA case, and it remains 
below IceCube observability in the BWCA case. The rate for neutralino annihilation in 
the sun or earth is given by 

T A = ^Ctanh 2 (i /r), (3.7) 



- 19 - 



where C is the capture rate, A is the total annihilation rate times relative velocity per 
volume, t & is the present age of the solar system and r = \j\fCA is the equilibration 
time. For small C, as shown in frame a), the equilibration time becomes large, so that 
~ 7}C 2 At 2 , and is hence sensitive to the neutralino annihilation cross section times 
relative velocity, unlike cases where the neutralino-nucleon scattering cross section is large. 
The muon flux jumps to observable levels at more negative values of Mi, but this is only 
because the Z\ suddenly becomes wino-like, so that the relic density becomes too low. 

In frames c), d), e) and f) we show the flux of photons, positrons, antiprotons and 
antideuterons, respectively. The results here are plotted as ratios of fluxes normalized to the 
mSUGRA point, in order to give results that are approximately halo-model independent. 
(We do show the above described expected experimental reach lines as obtained by using 
the Adiabatically Contracted N03 Halo model [50].) The rates for indirect detection via 
observation of halo annihilation remnants are typically low in the BWCA scenario, since 
the bino annihilation cross sections are s-wave suppressed. Observable results are indicated 
for 7 rays by GLAST, but this is due in part to the very favorable N03 halo distribution 
which is assumed. 



In Fig. 11, we show the same direct and indirect DM detection rates as in Fig. 10 



except this time versus Mi instead of M\. In frame a), the neutralino-nucleon scattering 
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Figure 10: Rates for direct and indirect detection of neutralino dark matter vs. Mi for mo = 
m i/2 = 300 GeV, with tan/3 = 10, Aq = 0, [i > 0. Frames c) -f) show the ratio of indirect detection 
rates compared to the mSUGRA model. In this plot, we adopt the adiabatically contracted N03 
distribution for halo dark matter. 
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rates do not have negative interference, and can remain at observable, although not en- 
hanced, levels. The rates for detection of BWCA DM at IceCube are relatively low, as are 
rates for anti- matter detection by Pamela in frames d) and e). The rates for 7 detection 
by GLAST in frame c) are similar to those from the mSUGRA case, while the rate for 
antideuteron detection by GAPS is just barely observable for BWCA DM in frame j). 

Overall, prospects for direct or indirect detection of BWCA dark matter are generally 
at or below levels expected in the mSUGRA model. For this reason, we do not present 
direct and indirect detection rates in the tuq vs. my 2 plane. We also mention that the 
situation is sharply different in the MWDM scenarios where the corresponding rates are 
generally larger than in the mSUGRA model. Thus, a detection of a signal in the XENON 
or Pamela experiments could serve to discriminate between these scenarios, especially if 
we already have some information of the SUSY spectrum from collider experiments. 



4. Neutralino radiative decay in the BWCA scenario 

The loop-induced radiative decay width for Z2 — > ^17 has been calculated in Ref. [80 



81]. A thorough numerical analysis p2, 83] has shown that the radiative decay can be 



large and even dominant in certain regions of MSSM parameter space. A necessary (but 
not sufficient) condition for this is that all tree level two body decay modes of Z2 be 
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Figure 11: Rates for direct and indirect detection of neutralino dark matter vs. M2 for mo = 
m i/2 = 300 GeV, with tan/3 = 10, Aq = 0, \i > 0. Frames c) -f) show the ratio of indirect detection 
rates compared to the mSUGRA model. In this plot, we adopt the adiabatically contracted N03 
distribution for halo dark matter. 
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kinematically forbidden. In this case, Z2 usually decays via Z2 — * Ziff, where / is a 
light SM fermion. However, if Z\ and Z2 are close in mass, the formally higher order two 
body radiative decay becomes competitive with the three body decays Z2 — > ffZ\. This 
is because T(Z2 — > -Z'it) oc (1 — /mg while T(Z2 — > Z\ff) oc (1 — m^/mg ) . 

As we have seen, in both the MWDM and BWCA scenarios — is small, so 
that we may expect that the branching fraction for radiative decays may be enhanced. 
Moreover, in both cases, the couplings of the neutralinos to the Z boson, which occur only 
via the higgsino components of the neutralino are strongly suppressed so that virtual Z 
boson exchange contribution to three body decay amplitudes is correspondingly suppressed. 
However, vector boson-gaugino loops essentially also decouple from the radiative decay in 
the BWCA case because the bino does not couple to these. These do not, however, decouple 
in the MWDM case since both the photino and the zino couple to W^W^ system. The 
branching fraction for the radiative decay is thus a result of a complicated interplay between 
the kinematic and dynamic effects discussed above. 



In Fig. 12, we show in upper frame a) the neutralino relic density, and in b) the 
BF(Z2 — * Zij) for the same parameters as in Fig. ffl, versus GUT scale gaugino mass Ml. 
The radiative branching fraction at this point in the mSUGRA model is just ~ 10 -5 . As 
Mi climbs to 490 GeV, in the MWDM case, the branching fraction has climbed to ~ 5%. 
When Mi varies to large negative values in the BWCA case, the branching fraction has 
climbed to 10%. In this case, we may expect that a considerable fraction of SUSY events 
at colliders to contain hard isolated photons via the decay of Z2 that is produced either 
directly, or via cascade decays of heavier sparticles. In the two lower frames, we show the 
same figures except for large mo = 1 TeV. In this case, the sfermion loops mediating the 
Z2 — ► Z\"f decay become suppressed, and the branching fraction is much smaller, reaching 
~ 8% in the MWDM case, and just 1% in the BWCA case. 

A similar situation is illustrated in Fig. [l^, where now we plot versus variable M2, 
while keeping Mi = = 300 GeV. In the upper frames for tuq = 300 GeV, we see 

that while BF(Z 2 Z\i) reaches 5% in the MWDM case, it reaches 25% in the BWCA 
scenario. In the lower frames, we see that for mo = 1 TeV, the branching fraction reaches 
8% for the case of MWDM and 0.8% for BWCA dark matter. We note that in both figures 
B(Z2 — ► Zij) attains a higher value at its peak when M1/M2 is negative because it is in 
this case that the sparticle masses get really close (see the level crossings in Figs. || and |||). 



In Fig. 14, we plot contours of BF{Z2 — > #17) in the tuq vs. m^M plane for Aq = 0, 
tan/3 = 10, n > and mt = 178 GeV. At every point in the plane of frame a), we have 
adjusted Mi to a negative value chosen so that tig h 2 = 0.11 in the BWCA scenario. We see 
that the BF(Z 2 -» Z17) exceeds 50% around m ~ 600-1000 GeV for m 1/2 ~ 1-1.4 TeV. 
The branching fraction remains large at all m^j values, but diminishes for mo ~ 500—1000 
GeV. In Fig. 06), we plot the same contours except that at every point in the plane, we 
have instead adjusted M2 to negative values so that ^j^ 2 = 0.11. In this case, we see 
again that the branching fraction can become larger than 0.3 for low mo values. In frames 
c) and d) once again we show contours of B(Z 2 — * Zyy) but for the case of MWDM, where 
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m =1000 GeV, m 1/2 =300 GeV, tanp=10, A Q =0, \i>0, m=178 GeV 




M 1 (GeV) M 1 (GeV) 



Figure 12: Upper: A plot of a) relic density QcDNih 2 and b) BF(Z2 — ► ^17) as a function of M\ 
for too = 300 GeV, to 1/2 = 300 GeV, A = 0, tan/3 = 10, /x > and m f = 178 GeV. Lower: Same 
plot for too = 1000 GeV. 

Mi and M2 are dialed to positive values. While the radiative branching fractions never 
reach much beyond the 10% level for MWDM, they maintain a significant rate out to large 
values of mo- This is because in the MWDM case the radiative loop decays are dominated 
by W-chargino exchange, whereas in the BWCA case the radiative loops are dominated by 
fermion-sfermion exchange. 

5. BWCA dark matter at colliders 
5.1 Z2 — Z\ mass gap in BWCA scenario 

An important question is whether collider experiments would be able to distinguish the 
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Figure 13: Upper: A plot of a) relic density VLcDKih 2 and b) BFiZ^ — > Z17) as a function of M2 
for mo = 300 GeV, to 1/2 = 300 GeV, A = 0, tan/3 = 10, /x > and m t = 178 GeV. Lower: Same 
plot for too = 1000 GeV. 



case of BWCA dark matter from other forms of neutralino DM such as MHDM as occur 
in the mSUGRA model or MWDM. We have seen from the plots of sparticle mass spectra 
that the squark and gluino masses vary only slightly with changing M\ or M2. However, 
the chargino and neutralino masses change considerably, and in fact rather small mass gaps 
m w 1 ~ m Zi anc ^ m z 2 ~ m Zi are ^ n S enera l expected in both BWCA and MWDM scenarios, 
as compared to the case of models that incorporate gaugino mass unification close to the 
GUT scale. 



In Fig. 15, we show contours of the mass gap m^ 2 — in the tuq vs. my 2 plane 
for A = 0, tan /3 = 10 and fi > for a) the mSUGRA model, b) the case of BWCA 
DM where —M\ is raised at every point until $7^ h 2 — ► 0.11 and c) the case of BWCA 
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BWCA: M 1 * m v2 , tan$=10, A =0, \i >0, m t =178 GeV 



BWCA: M 2 * m 1/2 , tan$=10, A =0, n >0, m, =178 GeV 




MWDM: M t * m 1/2 , tan$=10, A =0, \i >0, m, =178 GeV 
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Figure 14: Contours of BF{Zi — > Zij) in the mo vs. mi/? plane for tan/3 = 10, Ao — 0, fi > 0, 
m t = 178 GeV. In a)., we have adjusted Mi everywhere to negative values so that fig /i 2 = 0.11 in 
the BWCA scenario. In &)., we have adjusted Mi everywhere to negative values so that fij^/i 2 = 
0.11 in the BWCA scenario. In c). (d).), we have adjusted Mi and (M2) everywhere to positive 

values so that fl? h 2 = 0.11 in the MWDM scenario. 

^1 



DM where — M2 is lowered until ^^^h 2 — > 0.11. In the case of the mSUGRA model, most 
of the parameter space has mg — > 90 GeV, which means that Z 2 — ► Z\Z® decay 
is allowed. When this decay is allowed, its branching fraction is always large, unless it 
competes with other two-body decays such as Zi — > Z\h or Z2 — * ff or // (where / is 
a SM fermion). In the case of BWCA DM in frames b) and c), we see that (aside from 
the left-most portion of frame b), which is not a region of BWCA), the mass gap is much 
smaller, so that two-body tree level decays of Z2 and W\ are closed and three-body decays 
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mSUGRA: tanf,=10, A =0, p. >0, m, =178 GeV 
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Figure 15: Contours of — mass gap in the too v s- m i/2 plane for tan/3 = 10, j4o = 0, 
fi > and a) mSUGRA model, 6) -Mi > m 1/2 BWCA and c) -M 2 < m 1/2 BWCA. 



are dominant. If i^'s are produced at large rates either directly or via gluino or squark 
cascade decays^], it should be possible to identify opposite sign/ same flavor dilepton 
pairs from their decays, to reconstruct their invariant mass, and extract the upper edge of 
the invariant mass distribution [K c 



Finally, we note one curious feature of Fig. |l5|c that was referred to in the caption of 
Fig. |(| Within the blue shaded LEP2 excluded region there appear three allowed strips. 
The lower horizontal strip at m^ 2 ~ 100 GeV corresponds to the neutralino Z-annihilation 
funnel, where M2 does not have to be dialed to low values, since Z annihilation already 
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reduces the relic density. However, in this case where M2 is negative, terms in the two-loop 
gaugino mass RGE conspire to yield a weak scale \M%\ value which is somewhat larger than 
its mSUGRA counterpart. This means that while 2m^ ~ Mz, at the same time tti^ is 
just above 103.5 GeV, and is thus LEP2 allowed. A similar situation occurs for the longer 
allowed strip at ~ 150 GeV, except in this case it is the h resonance which reduces 
the relic density. Finally, at very low rriQ ~ 80 GeV, bulk annihilation via light sleptons 
reduces the relic density, while the negative M2 pushed above the LEP2 limit. The 
corresponding strips do not appear in frame b) because here, we have only scanned M\ 
values above mi/2, as required to get agreement with WMAP: as we can see from Fig. [j] 
that the Z and h resonances appear for \M\\ < m^. 

5.2 Fermilab Tevatron 

In the mSUGRA model, the best reach for SUSY at the Fermilab Tevatron occurs in 
the clean trilepton channel [CT, |88| . We examined the clean trilepton signal rate for case 
study point BWCA3. For this point, the total SUSY particle production cross section was 
~ 440 fb for y/s = 2 TeV pp collisions. Using the soft trilepton cuts SC2 of Ref. |88| 
(three isolated leptons withp T (4,4,4) > 11,7,5 GeV, |r/(£ lj2 / 3 )| < 1,2, EJp iss > 25 GeV, 
m(i£) > 20 GeV, mj-(f,i?y lss ) < 60 GeV plus a Z mass veto), we find a surviving signal 
cross section of only 0.035 fb, well below observability. The main problem here is the 
m(££) > 20 GeV cut. This cut is essential to remove background from virtual photons in 
qq' — > ivgl't' electroweak production. However, the cut also kills much of the SUSY signal 
in the BWCA case, since the it invariant mass is constrained to be < — , which 
is already quite small. 

Another possibility is to search for events at the Tevatron containing isolated high 
Et photons from radiative Z2 decays. We examined the isolated photon inclusive channel 
(E T (p() > 8 GeV with Y^cZe' 4 E T < O.IEt(t) and 177(7) | < 1.1), the photon plus lepton 
channel (p T (£) > 15 GeV with J2cZ° e A E T < 5 GeV and Ef iss > 20 GeV), and photon 
plus E^ iss > 20 GeV channels for the BWCA3 case study (see Table |). In the latter two 
channels, we also vetoed jets. The signal rates were 43, 2.7 and 13.5 fb respectively, while 
backgrounds from Wj (Z7) production were 6911 (1803), 2500 (6.2) and 394 (253) fb, 
respectively. In light of these results, it appears very difficult to use the isolated photon 
channels in the BWCA scenario to identify SUSY at the Tevatron. 

5.3 CERN LHC 

If the i?-parity conserving MSSM is a good description of nature at the weak scale, then 
multi-jet plus multi-lepton plus E™ 1SS events should occur at large rates at the CERN LHC, 
provided that rrig ~ 2 — 3 TeV[jS9|. The LHC reach for SUSY in the mSUGRA model has 
been calculated in Ref. ^TJ. The ultimate mSUGRA reach results, coming from the 
£;miss_|_ j e £ g channel should also apply qualitatively to the BWCA case, since the values 
of rrig and rriq change little in going from mSUGRA to BWCA, and the E™ ss + jets reach 
mainly depend on these masses. 

The reach in other channels such as multileptons plus jets and isolated photons plus jets 
may change substantially in the BWCA scenario. The reach of the LHC in the (mo, mi/2) 
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plane of the mSUGRA model was recently re-assessed in Ref. |9l|] . The search strategy was 
based on the detection of gluino and squark cascade decay products, namely multiple high 
transverse momentum jets and/or leptons and/or photons plus large missing transverse 
energy. Here, we use Isajet 7.72' [^] for the simulation of signal and background events at 
the LHC. The event and detector simulation was performed along the lines established in 
Ref. [31 1, where details on cuts and detector resolution along with our definitions of jets 
and isolated leptons and photons may be found. 

We plot the reach of the LHC 



BWCA: M 2 * m 1/2 , tan$=10, A =0, \i >0, m t =178 GeV 
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Figure 16: Reach contours of the CERN LHC for 100 
fb _1 of integrated luminosity for various signal topologies 
in the mo vs. m±/2 plane for tan (3 — 10, Aq = 0, /i > and 
for —M<z < m 1 / 2 such that f^^i 2 = ^-H a t every point in 
the plane. 



in Fig. 16 for the BWCA case us- 
ing the procedure described in |)]J . 
All events had to pass the pre-cuts, 
which impose the requirement that 
Ef iss > 200 GeV and there are at 
least 2 jets with > 40 GeV. 
We then optimized the cuts using 
the strategy in Ref. [91] - gener- 
ally speaking, harder jet Et and 
E™ ss cuts apply for heavier spar- 
ticles, and softer cuts apply for lighter 
sparticles. The events are divided 
into several classes, characterized 
by the number of leptons or the 
presence of an isolated photon in 
the final state. The 5a discovery 
reach for 100 fb^ 1 of integrated lu- 
minosity is shown for the various 
channels. The ultimate reach in 
the E™ ss + jets channel ranges from 
m i/2 ~ 1-4 TeV at low mo to m^ ~ 
1 TeV at mo = 2 TeV. These re- 
sults are similar to those obtained 



in the mSUGRA model. The reach contour in the isolated photon plus jets plus S™ 1SS 
channel, however, has greatly increased in the BWCA scenario compared to the mSUGRA 
case. In mSUGRA, the reach contour varies from m-1/2 : 850 —* 500 GeV as mo : .1 — > 2 
TeV. In the BWCA case, where Z2 — * Zyy at a large rate, the photon reach contour reaches 
a maximum of m 1/2 = 1.2 TeV for m ~ 500 GeV. Thus, at CERN LHC, the BWCA sce- 
nario will be signalled by multijet plus isolated multilepton plus E™ 1SS events, but with a 
large content of hard isolated photons as well, at least for the case where mo ~ 1 TeV. 8 



8 Gauge mediated SUSY breaking (GMSB) models in which the next-to- lightest SUSY particle (NLSP) 
decays to a photon and a gravitino also have photons in SUSY events. It is unlikely that these models 
will be confused with the BWCA scenario because not only is the sparticle mass spectrum quite different 
but the photons in the GMSB scenario would typically have much larger energy because the gravitino is 
essentially massless. Moreover, unless the NLSP has many decay modes, the multiplicity of photons in the 
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For larger values of m,Q > 1 TeV, the Z 2 — * Zij branching fraction drops, and the reach 
projections become similar to the case of mSUGRA. We have also checked that the reach 
using the 7 + n lepton channel is smaller than the reach via the corresponding n lepton 
channel, primarily because the signal becomes too small to pass our 10 event/100 fb _1 
requirement. 

For SUSY searches at the CERN LHC, Hinchliffe et al. have pointed outQ that an 
approximate value of rriq or rrig can be gained by extracting the maximum in the M e ff 
distribution, where M eff = E™ iss + E T (jet 1) + E T (jet 2) + E T (jet 3) + E T (jet 4). Their 
analysis will carry over to the BWCA scenario, as well as in models with gaugino mass 
unification, so that the approximate mass scale of strongly interacting sparticles will be 
known soon after a supersymmetry signal has been established. 

In mSUGRA, a dilepton mass edge should be visible in SUSY signal events only if 
m 1/2 ~ 250 GeV or if Z 2 -> It, it decays are allowed. In the case of BWCA DM, as 
with MWDM, the dilepton mass edge should be visible over almost all parameter space. 
We illustrate the situation for four case studies listed in Table |2| 9 The first case, labeled 
mSUGRA, has ttiq = ro^ = 300 GeV, with Aq = 0, tan/? = 10 and [i > 0. In this 
case, gg, gq and qq production occurs with a combined cross section of about 12 pb, while 
the total SUSY cross section is around 13.4 pb (the additional 1.4 pb comes mainly from 
-ino pair production and -ino-squark or -ino-gluino associated production). The case of 
BWCA1, with Mi = —480 GeV, has similar rates of sparticle pair production. The case of 
BWCA2, with lighter chargino and neutralino masses, has a total production cross section 
of 19.2 pb, wherein strongly interacting sparticles are pair produced at similar rates as in 
mSUGRA or BWCA1, but -ino pairs are produced at a much larger rate ~ 6.1 pb. We also 
show the case of BWCA3, which is similar to that of BWCA2 except that [i < 0, which 
gives a better fit to (g — 2) M measurements. 

We have generated 50K LHC SUSY events for each of these cases using Isajet 7.72', and 
passed them through a toy detector simulation as described above. Since gluino and squark 
masses of the three case studies are similar to those of LHC point 5 of the study of Hinchliffe 



et al.\86\, we adopt the same overall signal selection cuts which efficiently select the SUSY 
signal while essentially eliminating SM backgrounds: E™ ss > max(100 GeV, 0.2M e ff), at 
least four jets with Et > 50 GeV, where the hardest jet has Et > 100 GeV, transverse 
sphericity S T > 0.2 and M eff > 800 GeV. 

In these events, we require at least two isolated leptons, and then plot the invariant 
mass of all same flavor /opposite sign dileptons. The results are shown in Fig. [l7]. In 
the case of the mSUGRA model, frame a), there is a sharp peak at m(£ + £~) ~ Mz, 
which comes from Z 2 — > Z\Z° decays where Z 2 is produced in the gluino and squark 



GMSB case would be much larger since every SUSY event would contain two photons. 

9 In this study, a toy detector simulation is employed with calorimeter cell size A77 x = 0.05 x 0.05 
and —5 < r\ < 5. The hadronic energy resolution is taken to be 80%/y/E for \rj\ < 2.6 and 100%/y/E 
for \r/\ > 2.6. The electromagnetic energy resolution is assumed to be 3%/y/E. We use a UAl-like jet 
finding algorithm with jet cone size R = 0.5 and p 1 ^ > 25 GeV. We also require that \r/i\ < 2.5 and 
< 3. Leptons (es or fjs) have to also satisfy py pton > 10 GeV. Leptons are considered isolated if the 
visible activity within the cone AR < 0.3 is Ei?y ells < 2 GeV. The strict isolation criterion helps reduce 
multi-lepton background from heavy quark (especially ti) production. 
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parameter 


mSUGRA 


BWCA1 


BWCA2 


BWCA3 


Mi 


300 


-480 


300 


300 


M 2 


300 


300 


-156 


-170 


A* 


409.2 


401.3 


402.3 


-401.6 


m § 


732.1 


733.4 


736.3 


736.8 


mix 


713.9 


715.3 


701.9 


703.5 




523.4 


535.0 


554.8 


566.4 




650.0 


651.8 


645.1 


646.1 




364.7 


371.0 


00 A n 

324.5 


327.7 




Q99 e 
0ZZ.0 


ooZ. ( 


399 R 


399 K 
OZZ.D 


m w 2 


439 Q 


49fi n 


41 Q fi 


491 3 


m m 


993 Q 


993 /I 
ZZo.4 


1 38 tZ 

loo.O 


1/11 7 
141. f 


m z 4 


433.7 


425.0 


/(ICO 

415.2 


/i 1 n T 

419.7 


m~ 

Zz 


414.8 


409.4 


414.0 


410.0 


m z 2 


223.7 


222.7 


138.6 


141.4 


m z, 


117.0 


200.0 


116.8 


118.8 


m A 


538.7 


537.1 


508.4 


508.5 


m H + 


548.0 


546.4 


517.9 


518.0 


m h 


115.7 


115.3 


114.0 


112.7 




1.1 


0.11 


0.10 


0.12 


BF(b -»■ 57) 


3.2 x 10~ 4 


3.3 x 10~ 4 


3.7 x 10~ 4 


4.6 x 10" 




12.1 x 10~ 10 


9.9 x 10- 10 


-13.8 x 10~ 10 


13.1 x 10" 


a sc (Zip) 


1.5 x 10~ 9 pb 


4.7 x 10~ 12 pb 


1.2 x 10~ 9 pb 


8.0 x lO" 11 


BF(Z 2 -► Z17) 


2.7 x 10~ 5 


0.096 


0.25 


0.25 



Table 2: Masses and parameters in GeV units for mSUGRA and three BWCA scenarios. In each 
case, m = m 1/2 = 300 GeV, A = 0, tan/3 = 10 and m t = 178 GeV. The third BWCA case has 
fi < 0, while the first two have fi > 0. 



cascade decays. In the case of BWCAl in frame b), we again see a Z° peak, although 
here the Z°s arise from Z3, Z4 and W2 decays. We also see the continuum distribution 
in m(l + l~) < — = 22.7 GeV. The cross section plotted here is ~ 0.035 pb, 
which would correspond to 3.5K events in 100 fb _1 of integrated luminosity (the sample 
shown in the figure contains just 135 events). In frame c)— with a cross section of ~ 0.07 
pb (but just 187 actual entries)- we see again the Z° peak, but also we see again the 
m(£ + £~) < m^ 2 —m^ = 21.8 GeV continuum. In both these BWCA cases, the — m^^ 
mass edge should be easily measurable. It should also be obvious that it is inconsistent 
with models based on gaugino mass unification, in that the projected ratios M\ : M2 : M3 
will not be in the order 1 :~ 2 :~ 7 as in mSUGRA. Although the Z2 — Z\ mass edge 
will be directly measurable, the absolute neutralino and chargino masses will be difficult 
to extract at the LHC. In frame d), we show the spectrum from BWCA3, which is similar 
to the case of BWCA2. 
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Figure 17: Distribution of same flavor/opposite sign dileptons from SUSY events at the CERN 
LHC from a) mSUGRA, b) BWCA1, c) BWCA2 and d) BWCA3 cases as in Table § 



5.4 Linear e + e collider 

The reach of the CERN LHC for supersymmetric matter is determined mainly by mq 
and rrig, which depend on mo and m\i2- In contrast, the reach of the ILC for SUSY is 
largely determined by whether or not the reactions e + e~ — > W^W^ or e + e~ — > £ + £~ 
are kinematically accessible [92]. For instance, chargino pair production is expected to be 
visible if y/s > 2m^ . The value of Tn^ depends mainly on Mi and /x. Thus, in the BWCA 
case where M<i = m\ii but M\ is variable, the reach of the ILC in the mo vs. m\i% plane 
will be similar to the case of the mSUGRA model. However, in the BWCA case where 
Mi = mj; 2 , with variable M2, the reach of the ILC will be enhanced compared to the 
mSUGRA case, since IM2I is typically much smaller for a given set of mo and m\/2 values. 
The situation is illustrated in Fig. [l8| where we show the ultimate reach of the LHC and the 
ILC in the mo vs. plane for tan/3 = 10, Aq = 0, \i > and m t = 178 GeV. We have 
dialed —M2 at every point to give Q^h 2 = 0.11, in accord with the WMAP observation. 
We have assumed 100 fb _1 of integrated luminosity for both LHC and ILC. The reach 
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of ILC with y/a = 500 GeV (denoted by ILC 500) extends to m l/2 ~ 500 GeV, while 
the corresponding reach in the mSUGRA model with gaugino mass unification extends to 
m 1/2 ~ 320 GeV@. The reach of ILC with y/s = 1000 GeV extends to m 1/2 



1 TeV, 

compared with the mSUGRA value of mu2 ~ 600 GeV. In fact, we see that for mo ~ 1.9 
TeV, the ILC 1000 reach begins to exceed that of the LHC. In this region, ~ 500 GeV, 

2675 GeV. 



while nig 



2300 GeV and m~ q - 



BWCA: M 2 * m 1/2 , tan$=10, A =0, n >0, m t =178 GeV 
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Figure 18: Reach of the ILC and the LHC for SUSY 
in the BWCA DM scenario where \M2\ is lowered until 



fly h 2 = 0.11 at every point in the mo vs. 



l l/2 



plane. 



We assume 100 fb 1 of integrated luminosity for each 
collider, and show the ILC reach for yfs = 500 GeV 
and 1000 GeV. 

cation is realized in nature. 



At a y/s = 500 GeV ILC, the new 
physics reactions for the four case stud- 
ies shown in Table Q would include 
Zh, Wi"Wi, Z1Z2 and Z2Z2 produc- 
tion. It was shown in Ref. J92[| that 
even in the case of a small W\ — Z\ 
mass gap chargino pair production events 
could still be identified above SM back- 
grounds using cuts specially designed 
to pick out low visible energy acollinear 
signal events over backgrounds from 
e + e~ and 77 processes. The chargino 
and neutralino masses can be inferred 
from the resultant dijet distribution in 
WfWi -> {Iu i Z 1 ) + (qqZi) events^, 
pi) , p^ ]. These measurements should 
allow the absolute mass scale of the 
sparticles to be pinned down, and will 
complement the Z2 — Z1 mass gap mea- 
surement from the CERN LHC. The 

m wv m z 1 and 
nig measurements will point 

to whether or not gaugino mass unifi- 



combination of m^ 2 



m z 2 



The dependence of the cross sections for W-j W± , Z1Z2 and Z2Z2 production on the 
longitudinal polarization of the electron beam provides an additional tool at the ILC. This 
dependence is illustrated in Fig. 19 for the mSUGRA model in Table ^ in frame a), for case 
BWCAl in frame b) and for case MWDM1, which is point 1 of Ref. [24]. In all three cases, 
the chargino and the second lightest neutralino have dominant wino components so that 
the magnitudes and the polarization depence of cr(VUiVUi) and (7(^2^2) are qualitatively 
similar to one another and to the corresponding dependence in the mSUGRA modelp4|]. 
A minor difference is that, for Pl = —1, (7(^2^2) does not fall to quite as small values for 
the MWDM1 case because its hypercharge gaugino component always remains significant, 
as we have already discussed. The polarization-dependence and/or the magnitude of the 
e + e _ — ► Z1Z2 process are, however, quite different in the three cases. In the limit that 
Z\ is a bino and Z 2 is the wino (a good approximation in the mSUGRA case and an 
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Figure 19: Cross sections for e + e — > W^W Y , Z\Zi and Z2Z2 processes versus electron beam 
polarization Pl{c~) for an electron-positron linear collider operating at s/s — 500 GeV for a) 
mSUGRA, 6) BWCAl and c) MWDM1 from Ref. @. The positron beam is taken to be unpolar- 
ized. 



even better approximation in the BWCA case), i-channel selectron exchange dominates the 
amplitude. 10 Recall, however, that for Pt,(e~) = 1 {Pi{e~) = —1) only e/, (e^) exchange is 
possible. Since the chargino is always wino-like, its couplings to e# are strongly suppressed, 
accounting for the behaviour of the cross sections in the first two frames. For the MWDM1 
case in frame c), we would naively expect that Z\ would be photino-like and Z2 would be 
zino-like. It turns out, however, that (for the specific MWDM1 parameters) the difference 
Afi(weak) — M2(weak) is sufficiently large, and results in a flip of the relative sign between 
the gaugino components of both Z\ and Z2, while roughly preserving the magnitude. This 
strongly suppresses the eeiZi coupling, and hence, the exchange amplitude, resulting 
in the relatively flat and even decreasing polarization-dependence of the cross section. 11 

A very striking feature of the figure is the very large cross section for Z1Z2 production 
for a left polarized electron beam for the BWCAl case in frame b), as compared with 
the mSUGRA case. In both cases, since Z\ and Z2 couple to (a not-so-heavy) e/, mainly 
via their large hypercharge and SU(2) gaugino components, respectively, ei exchange 
completely dominates Z1Z2 production at Pl = 1- Moreover, the magnitudes of the 
couplings, as well as selectron masses, are very comparable in the two scenarios. The 
reason for the difference in the cross sections lies in the relative sign between the Z\ and 
Z2 eigenvalues (not physical masses) of the neutralino mass matrix. Since M1/M2 < 0, we 
expect this sign to flip in the BWCA case as compared with the mSUGRA or the MWDM 
cases, where M1/M2 is positive. This is relevant because when we square the t-channel 
amplitude and sum over the neutralino spins, there is one term that is proportional to the 
product m z 1 m z 2 ( see e '9~ (^) °^ l^l> wnere the last term that includes the factor 



(— is the one we are referring to). This term always flips sign between any scenario 



10 Recall that Z couples to neutralinos only via their suppressed higgsino components. 
11 When both and cr exchanges are dynamically suppressed, the effect of the usually small Z exchange 
contribution (which leads to a more or less flat polarization dependence) may also be significant. 
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with positive M1/M2 vs. negative M1/M2. As a result, what was a significant cancellation 
in the mSUGRA case in frame a) of Fig. |l^ turns into a sum in frame b), accounting for 
the factor ~ 2.5 increase in a(Z\Z2) at Pl = 1. 

We examine this potential enhancement of the e + e~ — * Z\Z<i production cross section 



at a \fs = 0.5 TeV linear collider in Fig. 20. In frame a), we show contours of a{e + e 



Z1Z2) in fb in the mo vs. rrii/2 plane for Aq = 0, tan/3 = 10 and \i > 0, in the case of the 
mSUGRA model. The light (yellow) shaded region is where either chargino pair production 
or selectron pair production is kinematically accessible at a = 0.5 TeV ILC. We see that 
the Z1Z2 cross section only exceeds the 100 fb level at the very lowest values of mo and 
mi/2, in the lower left corner. The cross section at the ~ 10 fb level gives some additional 
reach of an ILC for SUSY beyond the range where chargino pair production or slepton pair 
production are possible pi], |92| . In frame b), we show the same situation for the BWCA 
scenario where M\ has been set to negative values everywhere in the mo vs. my 2 plane 
such that the WMAP value QgjT? = 0.11 is fulfilled. In this case, the Z1Z2 production 
cross section is in general increased everywhere, and exceeds 100 fb in a much larger region. 
Of course, the region of reaction kinematic accessibility has decreased somewhat due to the 
increase in m^ for a given mo and value, so that contours of accessibility reach only 
up to m.1/2 ~ 350 GeV (as opposed to 420 GeV in the mSUGRA case). Also, the cross 
section falls off at large values of mo since the amplitude is suppressed with increasing 
selectron mass. In frame c), we show the same cross section in the mo vs. my 2 plane, 
except here M2 is taken negative, and is decreased at each point in absolute value until 
the WMAP value is obtained. In this case, the kinematically accessible region increases 
since for a given mo and m^ 2 value m^ 2 is lowered compared to the mSUGRA case, while 

stays fixed. Again, a rather large region appears with a > 10 — 100 fb, so that Z1Z2 
production should be robust over much of parameter space in the BWCA scenario as long 
as Z1Z2 production is kinematically accessible and mo is not too large. An unexpectedly 
large value of a{Z\Z2) / (t{WiW\) may, therefore, be an indication that M1/M2 < and 
that the selectrons are not too heavy. 

The rather large e + e _ — ► Z1Z2 cross section at the ILC potentially leads to another 
distinguishing feature of SUSY events at the ILC in the BWCA model: at least for mo ~ 1 
TeV, we should expect a large rate for events with one or more isolated photons due to the 
enhanced branching fraction for Z2 — > Zij decays. Thus, from Z1Z2 production, we may 
expect 7+ ft events, while from Z2Z2 production, we may expect 77+ ft, 7 + jets+ ft 
events, 7 + 11+ ft events and 7+ ft. Since the 7 arises from a two body decay, the 
endpoints of the -£'(7) distribution will be functions of y/s, and ■ A measurement 
of the endpoints of this distribution will thus allow an independent measurement of the 
two lighter neutralino masses. The E{pf) distribution is shown in Fig. ^ for BWCA2 at a 
sfs = 0.5 TeV ILC. 

6. Conclusions 

In this paper, we have considered the phenomenological implications of neutralino dark 
matter in the bino-wino co-annihilation (BWCA) scenario, and compared these to the 
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mSUGRA: tan$=10, A =0, n >0, m, =178 GeV 
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Figure 20: Contours of er(e + e — > Z1Z2) in fb for a ^Js — 0.5 TeV ILC with unpolarized beams 
in the too vs. to^/2 plane for tan/? = 10, Aq = 0, fi > and a) mSUGRA model, b) —M± > m 1 / 2 
BWCA and c) —M 2 < m l / 2 BWCA. The yellow shaded regions show where W^W^ and/or 
should also be kinematically accessible to a y/s = 0.5 TeV ILC. 

case of mixed wino dark matter which has a qualitatively similar spectrum. The BWCA 
scenario arises in models with non-universal gaugino masses where, at the weak scale, 
I Mi I ~ I-M2I, but where the two gauginos masses have opposite signs. The same sign 
case gives rise to mixed wino dark matter, while in the BWCA case, there is little mixing 
so that the Z\ remains nearly a pure bino, while Z2 remains nearly a pure wino. The 
scenario can be brought into accord with the WMAP ^^/i 2 ~ 0.11 constraint by arranging 
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for rnTTr ~ 

Wi z 2 

annihilation mechanism in the early universe. 



~ Tag , so that bino-wino co-annihilation is the dominant neutralino 
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Figure 21: Distribution in E(j) from e + e — > Z2Z1 
with Z-2, — > 7Z1 for point BWCA2 in e + e~ collisions at 
= 0.5 TeV. 



serve to discriminate between these scenarios. 



Since co-annihilation processes are 
the dominant mechanism for the an- 
nihilation of relic neutralinos from the 
Big Bang, it is expected that indirect 
dark matter detection rates (which de- 
pend on the Z\Z\ pair-annihilation cross 
section) will typically be quite low in 
the BWCA scenario. Direct detec- 
tion of neutralinos may be possible 
for BWCA DM at stage 3 detectors 
in some portions of parameter space. 
In contrast, for the MWDM case, neu- 
tralino annihilation cross sections are 
enhanced relative to mSUGRA so that 
it is more likely that an indirect sig- 
nal for DM will be seen. If we already 
have some infomation about the spar- 
ticle spectrum from the LHC, it may 
be that results from direct and indi- 
rect DM detection experiments may 



The small mass gap expected between m^ 2 and is a hallmark of both the BWCA 
and the MWDM scenarios, and leads to a variety of interesting consequences for collider 
experiments. In teh BWCA case, if mo ~ 1 TeV, then the radiative decay Z2 — > Zij is 
greatly enhanced, leading to the production of isolated photons at the Tevatron, LHC and 
ILC. While extraction of the isolated photon signal above SM background looks difficult 
at the Tevatron, they should yield observable signals at the CERN LHC. The LHC should 
also be able to extract the m^—m^ mass difference from gluino and squark cascade decay 
events which contain isolated opposite sign/same flavor dilepton pairs. By comparing the 
neutralino mass difference to measurements of the gluino mass, we should readily be able 
to infer that the weak scale gaugino masses are incompatible with expectation from models 
with a universal gaugino mass at the high scale. At the ILC, again small fn^ — and 
m Z2 ~ m z± mass § a P s should be measurable, as in the MWDM and MHDM scenarios. It 
should also be possible to extract the weak scale values of the gaugino mass parameters. 
However, in the case of BWCA, Z1Z2 production should occur at enhanced rates relative 
to mSUGRA and MWDM, with a production cross section which increases with P/,(e _ ). 
In constrast to MHDM, the Z3 and Z^ states will typically occur at much higher mass 
scales associated with large values for the \n\ parameter. 
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Appendix: Effective Theories in Quantum Mechanics 

Consider a situation where the total Hamiltonian can be split into two pieces, 

H = H + V, 

so that the spectrum of Hq is hierarchical, i.e., it consists of "low energy" states 
possibly interacting with one another via interactions included in Hq, and "high energy" 
states, \h a ) (again possibly interacting with themselves via interactions included in Hq). 
Without loss of generality, all interactions between the states {\k}} and the states {|/i a )} 
are encapsulated in V. In other words, (h a \Ho\li) = 0, and {k\V\lj) = {h a \V\hb) = 0. 

We seek a description of low energy states in terms of an effective Hamiltonian il e ff, 
that acts only on the low energy states {|^)}, assuming that we are at energies that are 
too low to excite the high energy states. Note that even if Hq is completely diagonal in 
the low energy sector, we would expect that scattering of low energy states (with energy 
E <C M a ) would occur via their interactions with the high energy sector, with an amplitude 
suppressed by E/M a . Here, we show how this comes about and obtain an expression for 
H e s- 

We define H e g by matching the low energy Green's functions of the full theory with 
those of the effective theory, 

fcl ff+V J *i> = ^\jT—\h)> if 1*1 « M a- (1) 
riQ + V — z ti eS — z 

Here, z is the complex argument of the Green's function. Expanding the left hand side of 
(|l|), we obtain 

{^Itt^ — l z i> = ^Itt-i \h) + {k\-^— (-v)j^— \W 

H eS - z Hq-z J H - z H - z 

Ho — z Ho — z Ho — z 

Clearly, if V = (no interactions between low and high energy states), H e g is just Ho 
restricted to the low energy subspace, i.e. H c $ = HoV\ ow , where V\ ow is the projector 
on to the low energy subspace. It is also clear that terms with an odd number of (— V) 
factors on the right hand side of (||) are all zero because V only connects states in the low 
energy sector with those in the high energy sector, whereas Ho only connects low energy 
(high energy) states with one another. The third term on the right hand side of (||) can be 
written as 



H n — z J ,„., \ H — z J ab \Hq — a, j 
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with the repeated indices i',a,b and j' all summed up. We now defined an induced effective 
potential V indnced by, 

" (V induced ) = (-VU (-^—) (-V) bf , (3) 

V ,% 3 \ H 0- z Jab 



in terms of which this term can be written as 

' 5 ' yiadutvd 



H -zJ w \ h'j' \H -z, rj 

In exactly the same manner, the term with four factors of (— V) will end up as one with 
two factors of ^_y induced ^ e £ c . so that we have a geometric series that can be summed to 
give 

{kl H cS -J lj) = < Zi| ff + y induced _^>> W 

showing that the interactions between the low and high sectors effectively induce an addi- 
tional "potential" in the low sector. 

Up to now, our considerations (though formal) have been "exact" in the sense that 
we have not made any low energy approximation. This shows up in the fact that the 
"potential" y induced defined in (0) depends on the "energy" z. If \z\ is small compared with 
the high energy scale associated with the spectrum of Hq, we can ignore it in the evaluation 
of the matrix elements of (Hq — between high energy states, and approximate y induced 
by, 

(v*— % = - £ [v ia (^) ab v bJ + v ia o(K)v bj + 



a.b 



0/ ab 



where, in the last step, we have made the "low energy approximation" and obtained what 
is a conventional potential (independent of z). Assuming that the matrix elements V{ a have 
magnitudes corresponding to the low energy scale, we see that the (^jj^J term supresses 
the low energy matrix elements of y mduced by \z\/M a as expected. 
This expression is particularly useful in two cases. 

1. If HqV\ ow is diagonal, then all interactions arise only from y mduced an d this analysis 
is essential to obtain any scattering in the low energy sector. 

2. If the low energy sector has an approximate symmetry that is violated only by its 
interactions with the high mass sector, or even just by interactions solely within the 
high mass sector, we can use \/ mduced to study these symmetry violations. 

In the analysis up to now, we have retained just the leading correction in powers of M a . 
It is straightforward to retain the 0(1/M^) term in the expansion of the Green's function 
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(H c ff — z) 1 . The energy eigenvalues E in the low energy theory are given by those values 
of z where the corresponding Green's function develops a singularity, i.e., where 



H c r\ 



E\ 



Assuming, for simplicity, that Hq is diagonal in the high mass sector and retaining terms 
to 0(1/Mq), we find that these eigenvalues are given by the generalized matrix eigenvalue 
equation, 

j a a j a o. 

The positive operator on the right hand side serves as a metric in the low energy subspace. 

This formalism can be directly used to obtain the eigenvalues ( |2.9| ) of the neutralino 
mas matrix in the large |/x| limit. In this case, the low energy sector comprises of the 
neutral wino and the bino, and the generalized eigenvalue equation takes the form, 



M2 + a — atan#iy 
-a tan 8w Mi + a tan 2 9\y 

Ml 




- ¥fc tan 9 W 1 + tan 2 W 




(7) 



with a = sin 2/3 much smaller than M\ and M2 in the large limit. 
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